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We propose inelastic electron tunneling spectroscopy scanning tunneling microscopy (lETS-STM) 
as a means of exciting and observing intrinsic localized modes (breathers) in a macromolecule. As 
a demonstration, inelastic tunneling features of the density of states are calculated for a simple 
nonlinear elastic Morse chain. The general formalism we have developed for the lETS is applicable 
to other nonlinear extended objects, such as DNA on a substrate. 
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Dynamical localization phenomena have been a subject 
of intense theoretical research since discrete breathers (or 
intrinsic localized modes) were proven to be generic 
solutions of nonlinear spatially discrete systems. These 
solutions of the underlying equations of coupled nonlin- 
ear oscillators are characterized by being temporally pe- 
riodic and spatially localized y, Q . Two particularly 
important and physically appealing aspects of discrete 
breathers are 1) that they exist in translationally invari- 
ant systems and therefore these excitations are radically 
different from Anderson localization and other localiza- 
tion effects driven by external disorder or defects, and 
2) that they are very robust with respect to changes in 
the equation of motion. They have been discussed in 
such diverse systems as spin lattices , and dynamics of 
DNA HQ. 

The robust localization in the form of discrete 
breathers occurs because the discreteness of the system 
provides an effective cutoff for the wavelength of the lin- 
ear modes (e.g. phonons) that may exist in the system, 
and this can prohibit phonon resonances with all tem- 
poral harmonics of the discrete breather. This effectively 
eliminates decay mechanisms into extended linear modes. 
It is the nonlinearity of the systems that makes this possi- 
ble by admitting frequencies outside the phonon (linear) 
band. 

The field has been dominated by theoretical research, 
and general methods for direct experimental observa- 
tion of breathers are lacking, although experiments on 
localization of light propagation in weakly coupled op- 
tical waveguides 0, low-dimensional crystals and 
Josephson-junction networks 0, 0| have recently been 
reported. 

Here we propose the use of scanning tunneling mi- 
croscopy (STM) to both excite and observe breathers in 
an inelastic electron tunneling spectroscopy experiment 
(lETS). The lETS itself is a well estabhshed and power- 
ful tool that allows the measurement of the characteristic 
energies of local and extended modes. Examples of ap- 
plications of this technique, among many, include mea- 
surements of molecular stretching and vibrational modes 
in metal-insulator-metal tunnel junctions local 
lETS scanning tunneling microscopy (lETS-STM) mea- 



surements of the local vibrational mode of single atoms 
or molecules on a surface ! observation in tunneling of 
the collective magnetic resonance in the superconducting 
state of high-Tc materials O Ell; and recent lETS-STM 
of the spin flip of a magnetic atom in a magnetic fleld |l7l | . 
Here we give for the first time a theoretical basis for an 
lETS-STM observation of the breather mode in an ex- 
tended system. 
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FIG. 1: (Color) The proposed experimental realization: An 
extended molecule is placed on the metallic surface. An 
STM is used to inject /extract electrons from the molecule. 
Electron-vibration coupling produces excitations that are out- 
side of the phonon continuum. This localized excitation will 
re-form into a set of breather excitations with different en- 
ergies. The breathers will remain localized for a long period 
of time since they cannot decay into the phonon continuum. 
Subsequent lETS will reveal localized excitations. 

The lETS directly measures excitation energies. Elec- 
trons scatter off a local (or collective) mode. Due to 
the scattering, a contribution to the electron self-energy 
occurs above a corresponding threshold value of the fre- 
quency determined by the mode frequency. Thus, in a 
tunneling experiment, for a voltage bias exceeding the 
threshold, electrons can excite the mode. This additional 
scattering channel leads to a step in the density of states 
and in the tunneling conductance. Low temperatures 
are required to avoid thermal smearing of the step in the 
conductance. 

The setup we propose is to use STM in a two step 
process: First we propose to locally excite breathers by 
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injecting energetic electrons into the molecule on the sub- 
strate. The energy of the injected electrons is taken to be 
outside the phonon band and able to excite the breather; 
in the subsequent step we propose to measure the local 
density of states (LDOS) in the vicinity of the injection 
point. The first step is a major assumption of this work. 
There is no detailed microscopic theory for the energy 
transfer between electrons and breathers. We make the 
reasonable assumption that the resonance between the 
energy of the electron and that of the breather will al- 
low efficient energy transfer between electronic degrees of 
freedom and breathers. At the second stage, the LDOS at 
a finite bias but with much smaller drive currents will al- 
low measurement of the inelastic processes revealing the 
existence of the localized breather. Alternatively, one 
can view our work as addressing what will be the fea- 
tures in LDOS once the breathers are excited. The setup 
is illustrated in Fig. ^ 

In what follows we will use a simple model to demon- 
strate the possibility of lETS in macromolecules. We 
consider a translationally invariant long chain with vibra- 
tional sites n = — iV-f 1, . . . ,0, . . . , N with local molecular 
displacements x„ at each site, for simplicity taken only 
along the direction normal to the surface. The role of 
the surface in this consideration is to provide mechanical 
support for the chain and to provide conduction electrons 
whose DOS will ultimately be measured by STM. 

Specifically, we consider the following model Hamilto- 



H = 



n=N 

E 

n=-N+l 



(1) 



Here the first three terms describes the dynamics for the 
displacement field in an extended system, the fourth term 
is the Hamiltonian for uncoupled surface electrons with 
= ek — /i the energy dispersion measured with re- 
spect to the chemical potential, the last term Hint 



n=N 



is the interaction Hamiltonian that 



describes a local (Holstcin) coupling between the dis- 
placement Xn and local electronic density at a site n. We 
will ignore spin indices of the fcrmion operators c„ (cj^), 
since we are interested in a total electronic density and 
there is no explicit spin dependence in the Hamiltonian. 

We introduce the Matsubara Green's functions of the 
electrons on the substrate. The bare Green's function is 
given by: 

G°(k,zc^„) = , , (2) 



i2n) 



(3) 



where cun = (2n -I- l)7rT is the Matsubara frequency for 
electrons. In a two-dimensional (2D) system, assuming 



that the surface states of the electrons are interacting 
strongly with the displacements, we have G°(r,i(X'„) = 
—iTTNQSgn(ujn)Joikpr), with Nq being the density of 
states of the metal, Jo{x) the Bessel function of the zeroth 
order, and sgn(a;) the sign function. Employing pertur- 
bation theory up to second order (in the coupling con- 
stant A) in the presence of the dynamical disorder (dis- 
crete breathers) coupled to the electrons, the full Green's 
function is 

G(r, r';zw„) = G°{r,r';iujn) + J J dridr2G°{r,ri;iujn) 

(4) 



xE(ri, r2; iu;„)G°(r2, r'; iw„) 



with self-energy 
S(ri,r2;icj„) = 



X2 

-T 



-D°(ri,r2;il7„) 



xG°(ri,r2;i(w„ - 

(P'kdio np{^-\^ — 1 



(27r)2 iw„ - Ck - 
xS(ri,r2;cj). (5) 

Here _D"(ri, r2; ifim) is the Fourier transform 
of the correlation function _D°(ri,r2;ri — t-i) = 
(Tr [a:;(riri)a;(r2T2)]), fl„i = 2m7rT is the Matsubara fre- 
quency for the bosonic excitations, and the two distribu- 
tion functions are given by np^BiE) = l/[exp(i<^/T) ± 1], 
respectively. Of great importance in Eq. © is the 
generalized spectral function i3(ri,r2;w), which has not 
yet been specified and depends on the detailed model for 
the bosonic excitations. In this sense, Eq. Q together 
with Eq. (O is the most general formalism to describe 
the effects of the bosonic excitations on the electronic 
properties. 




FIG. 2: The displacement field in the breather excitation (a), 
and its Fourier transform b) loq = 0.7, A° - squares, A\ - 
circles, - small circles, A'^ - diamonds, Aj^ - pluses 

For the purpose of the present work, we need to identify 
the spectral function i3(ri, r2; i^) for the breather mode. 
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Since the dynamics for the breather in is described 
through a lattice model, we are actually studying the 
correlation function in discretized space coordinates. In 
real time and space, the retarded correlation function is 
given by: 

D°'''{n,m;t,t') ^ id{t ~t'){[x,-,{t),x,nit')]) , (6) 

where n,m £ molecule sites, i.e., r„ and r^, and Xn{t) is 
the displacement operator at the site n. We use standard 
quantized notations: 



a=l 



(7) 



with a being the eigenvalue index, corresponding to the 
energy of the mode uja = Qf^o with amplitude A" and 
ba being the mode annihilation operator. We determine 
the breather solutions as described in and obtain the 
amplitudes A" at site n along the molecule. The coeffi- 
cients A" describe the real space amplitude distribution 
of the breather mode with frequency tUa- This is a classi- 
cal solution that will not decay in time (we ignore weak 
damping to simplify the calculation). Then 



oo 

Q = l 



Its Fourier transform is obtained as: 

Aa* AO 



n ni 



Aa Aa* 



UJ + UJa + iS UJ — UJa + iS 



(8) 



(9) 



with infinitesimal 6, from which the spectral function fol- 
lows immediately as 




(b) 
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FIG. 3: (Color) The correction of the local density of states 
calculated at the center of the breather mode (a) and its 
derivative (b) as a function of electron energy for various tem- 
peratures T = 0.01 (Red), 0.05 (Green), and 0.1 (Blue). 



B{n,m;u;) = Y,A'^:A^[Siu; + u;^)-Siu;~LU^)] . (10) 



Here, without loss of generality, we have assumed A" to 
be real. The electronic self-energy in 2D is then found to 
be: 



E(r„, r„,; iw„) = X'^NoJoikprnm) J2a=i K.Ka 



Jde 



(11) 



These formulas are quite general. The only breather- 
specific features in Eqs. Ullfl and (|l()|l are that the 
breather amplitude factors Af^ are localized at the po- 
sition, say n = 0, where the breather is excited |l9l |. 

Once we arrive at the the electronic self energy due 
to the scattering of electrons off the breather mode that 
resides on the molecule, we are ready to calculate the 
change in the LDOS of electrons, which is given by: 



-IniG(r„r„S-HiO+) 



(12) 



where the retarded Green's function is obtained through 
the analytic continuation from Eq. Q 

G{ri,rf,E + iO+) ^ G"(0; E + iO+) 
+ En.m G"(r. - r„; S + zO+)E(r„, r„; S + zO+) 

xGO(r„ -r2;£; + iO+) . (13) 

The LDOS is proportional to the low temperature lo- 
cal tunneling conductance, which can be measured by 
STM [13. 

A little algebra leads to the correction to the LDOS: 



SN{Ti,E) 
No 



XJo{kFrin)Jo{kF )Jo{kFr„ii) 
x[nF{E + LUa) - nF{E - UJa) + nsiuJa) - nBi-LOa)] . 

(14) 

Equation (|14|l is the central result of this paper. It allows 
one to calculate the corrections to the LDOS at any point 
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FIG. 4: (Color) The spatial image of the derivative of the 
LDOS correction at energy E = ujq = 0.7 due to the electron 
scattering off the breather mode. The temperatures T = 0.01. 



on the surface and quantify the effects of the inelastic 
tunneling corrections due to coupling to a breather. 

For the numerical calculation, we take the parameter 
values in Eq. The atomic mass M = 1, the elas- 
tic constant k = 0.1, and a Morse potential V{x) = 
D[l-exp{-x^/^^)] with D = 0.5 and = 1, so that the 
bottom of the phonon continuum is located at Uc = 1. 
The lattice constant of the chain molecule is taken to be 
^0- Energy is measured in units of Uc and length in units 
of ^0 ■ We then choose the fundamental frequency of the 
breather mode loq = 0.7, and the dimensionless electron- 
mode coupling constant (AiVo)^ = 10~^. As shown in 
Fig. 121 the breather excitation is localized at some region 
in the molecule with a length scale of about four lattice 
constants along the chain direction, in contrast to the vi- 
bration mode of a single-atom molecule, which is located 
on a single site. Moreover, due to the anharmonicity of 
the Morse potential, the coefhcients for higher frequency 
Fourier modes are finite. All these unique features can 
be detected by the lETS-STM experiments. Figure H il- 
lustrates the correction to the LDOS (upper panel) due 
to the presence to the breather, and the effects on the 
derivative of the LDOS (lower panel). As expected, there 
are steps in the LDOS occurring at the energy E = ujq, 
2u}o, Swoi • ■ • ■ For example the step at E = 2ujq originates 
from the inherent anharmonic nature of the breather but 
has much weaker intensity than the first step. However, 
corresponding to these steps, the peaks in the derivative 
of the LDOS are easier to identify. These features are 
however rapidly smeared out as the temperature is in- 
creased. Therefore, to detect the breather excitation, 
the temperature must be much lower than the mode 



frequency. Experimentally, the temperature of several 
Kelvins and a typical value of ujq ~ lOOmeV are easily 
accessible. 

Figure 01 shows the spatial distribution of the deriva- 
tive of the LDOS correction at energy E = ojq. The 
strongest intensity occurs at the center of the localized 
breather core. However, one can also clearly see the spa- 
tial extension away from the breather (especially along 
the chain direction lying horizontally), which can be de- 
tected again by the lETS-STM. The ripples in the far 
field from the breather comes from the rapidly oscillat- 
ing factors Jo{kfr). All these features are observable by 
the lETS-STM by measuring the local tunneling conduc- 
tance dl/dV and its derivative (Pl/cPV. The strongest 
signature of the localized lattice excitation is in the vicin- 
ity of the breather. 

In summary, we have proposed a novel application of 
STM to perform lETS on extended molecules. Intrinsi- 
cally localized nonlinear excitations (breathers) can first 
be excited and subsequently imaged using lETS STM. 
We find that indeed breathers result in an extremely lo- 
calized lETS features and may be observed with STM 
at low temperatures. The formalism we have developed 
for the lETS of an extended objects is quite general and 
is applicable to other systems. One of the most inter- 
esting extended systems where this experiment could be 
performed is lETS on DNA. 
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